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o : 

O ■ We provide a theoretical study of the quantum adiabatic evolution algorithm with different 

n 

p ^. evolution paths proposed in |1[]. The algorithm is applied to a random binary optimization 

ly^ [ problem (a version of the 3-Satisfiability problem) where the n-bit cost function is sym- 

^ ' metric with respect to the permutation of individual bits. The evolution paths are produced, 

^ . using the generic control Hamiltonians H{t) that preserve the bit symmetry of the under- 

> ; 

^ , lying optimization problem. In the case where the ground state of H{0) coincides with 

in ■ 

I ' the totally-symmetric state of an n-qubit system the algorithm dynamics is completely de- 

o\ : 

O ■ scribed in terms of the motion of a spin-n/2. We show that different control Hamiltonians 

m ■ 

O I can be parameterized by a set of independent parameters that are expansion coefficients of 

H{t) in a certain universal set of operators. Only one of these operators can be responsible 
P5 ■ for avoiding the tunnelling in the spin-n/2 system during the quantum adiabatic algorithm. 

^ : 

^ . We show that it is possible to select a coefficient for this operator that guarantees a poly- 

^ I nomial complexity of the algorithm for all problem instances. We show that a successful 



^ ' evolution path of the algorithm always corresponds to the trajectory of a classical spin-n/2 

and provide a complete characterization of such paths. 

PACS numbers: 03.67.Lx,89.70.+c,75.45.+j 



I. INTRODUCTION. 

Recently a novel paradigm was suggested for the design of quantum algorithms for solving 
combinatorial search and optimization problems based on quantum adiabatic evolution In the 
quantum adiabatic evolution algorithm (QAA) a quantum state is closely following a ground state 
of a specially designed slowly time- varying control Hamiltonian H{t). At the beginning of the 
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algorithm the control Hamiltonian H{0) = Hb has a simple form with a known ground state that 
is easy to prepare, and at the final moment of time it coincides with the "problem" Hamiltonian 
Hp which ground state encodes the solution of the classical optimization problem in question 



ffp = ^E,|z)(z| (1) 

z 

|z) = \zi)i (g) 12:2)2 ® ■ ■ • ® \Zn)n- (2) 

Here E^, is a cost function defined on a set of 2" binary strings z = {zi, . . . , z„} z^ = 0, 1, each 
containing n bits. The summation in ([T]) is over the 2" states |z) forming the computational basis 
of a quantum computer with n qubits. State \zk)k of the k-th qubit is an eigenstate of the Pauli 
matrix with eigenvalue 1 — 2zk ±1. If at the end of the QAA the quantum state is sufficiently 
close to the ground state of Hp then the solution to the optimization problem can be retrieved by 
the measurement. ^ 

It has been shown recently |7] that the query complexity argument that lead to the exponential 
lower bound for the unstructured search ^ cannot be used to rule out the polynomial time solution 
of NP-complete Satisfiability problem by the quantum adiabatic evolution algorithm (QAA). 

A set of examples of the 3-Satisfiability problem has been recently constructed Q] to test 
analytically the power of QAA. In these examples the cost function depends on a bit-string z 
with n bits, z = {zi , Z2, ■ ■ ■ , Zn}, only via a Hamming weight of the string, = zi + Z2 + ■ ■ ■ + Zn, 
so that Ez = f (wz) where the function f{w) is in general non-monotonic and defines a particular 
instance of this "Hamming Weight Problem" (HWP). In [0,0] the original version of QAA 10] 
was applied to the HWP where the control Hamiltonian is a linear interpolation in time between 
the initial and final Hamiltonians. 

In this case, it was shown that the system can be trapped during the QAA in a local 

minimum of the cost function for a time that grows exponentially in the problem size n. It was 
also shown 1I4 ] that an exponential delay time in the quantum adiabatic algorithm can be inter- 
preted in terms of the quantum-mechanical tunnelling of an auxiliary large spin between the two 
intermediate states. 

The above example has a significance greater than just being a particular simplified case of a 
binary optimization problem with symmetrized cost. Indeed, one can a rgu e that it shows a generic 
mechanism for setting "locality traps" in the 3-Satisfiability problem [ 10]. But most importantly, 
this example demonstrates that exponential complexity of QAA can result from a collective phe- 
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nomenon in which transitions between the configurations with low-lying energies can only occur 
by simultaneous flipping of large clusters containing order-n bits. In spin glasses, there is typically 
an exponential number of such configurations, the so-called local ground states. A similar picture 
may be applicable to random Satisfiability problems 1 15]. In some cases, these transitions can be 
understood and described in terms of macroscopic quantum tunnelling. A tunnelling of magne- 
tization was observed in large-spin molecular nanomagnets lllll and in disordered ferromagnets 

Q. 

n 

The paper [1] suggests that large tunnelling barriers can be avoided in QAA by using multiple 
runs of QAA with realizations of the control Hamiltonians H{t) sampled from a random ensem- 
ble. This ensemble is chosen in a sufficiently simple and general form that does not depend on the 
specific instance of the optimization problem. Different Hamiltonians H{t) correspond to differ- 
ent paths of the unitary evolution that begin and end in the same initial and final states (modulus 
phase factors). The complexity of QAA with different paths for the HWP was tested numerically 
in jj] using an ensemble of random 8x8 matrices. The results indicate that the HWP may be 
solved in polynomial time with finite probability. 

In case when the random paths H{t) preserve the bit-permutation symmetry of the problem 
it is natural to describe the random ensemble of H{t) in terms of the dynamics of a spin-r2/2 
system. This approach allows for a general theoretical analysis of the algorithm. In the present 
paper, we perform this analysis for the random version of HWP (over-constrained 3-Satisfiability 
problem) by mapping the dynamics of QAA onto the motion of a quantum particle in a ID effective 
potential. This allows us to compute the statistical weight of the successful evolution paths in the 
ensemble and provide a complete characterization of such paths. 

II, QUANTUM ADIABATIC EVOLUTION ALGORITHM WITH DIFFERENT PATHS 

In a QAA with different paths 111, one specifies the time-dependent control Hamiltonian 

H{t) ^ H{t) 

H{t) = {l-r)HB + r(l -r)HE + r Hp, (3) 

r=^ e(0,i). 

where the control parameter r plays the role of dimensionless time. This Hamiltonian guides the 
quantum evolution of the state vector \^{t)) according to the Schrodinger equation ih d\'ijj{t))dt = 
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H{T)\tp(t)) from t = Otot = T, the run time of the algorithm. Hp is the "problem" Hamiltonian 
given in ([T]). Hb and He are 'driver" Hamiltonians designed to cause the transitions between the 
eigenstates of Hp. 

An initial state of the system |V'(0)) is prepared as a ground state of the initial Hamiltonian 
if(0) = Hb- It is typically constructed assuming no knowledge of the solution of the classical 
optimization problem and related ground state of Hp. In the simplest case 

n 

HB = -CY.al |^(0)) = 2-"/2j]|z), (4) 

j=l z 

where cr^ is a Pauli matrix for j-th qubit and C > is some scaling constant. The ground state of 
Hb has equal projections on any of the 2" basis states |z) ©. 

Consider instantaneous eigenstates \4>k{T)) of H{t) with corresponding eigenvalues Ek{T) ar- 
ranged in non-decreasing order at any value of r G (0, 1) 

H{r)\cPk{r)) = A,(r)|0fc(r)), A; = 0, 1, . . . , 2" - 1. (5) 

Provided the value of T is large enough and there is a finite gap for all t G (0, T) between the 
ground and exited state energies, AA(r) = Ai(r) — \q{t) > 0, quantum evolution is adiabatic 
and the state of the system stays close to an instantaneous ground state, |0o(^/^)) (^P to a 

phase factor). Because H{t) = Hp the final state |^/'(T)) is close to the ground state \(f)o{r = 1)) 
of the problem Hamiltonian. Therefore a measurement performed on the quantum computer at 
t = T will find one of the solutions of combinatorial optimization problem with large probability. 
Quantum transition away from the adiabatic ground state occurs most likely in the vicinity of the 
point T ^ Tc where the energy gap AA(r) reaches its minimum (avoided-crossing region). The 
probability of the transition is small provided that IitI 

T> ^lij^axAA"!, (6) 

where 

dH 

iJinax = max K0i(r)| — |0o(i"))|, 
rG(o,i) ar 

AAmin = min [Ai(r) - Ao(r)] , (7) 

Te(o,i) 

The r.h.s. in Eq. Q gives an upper bound estimate for the required runtime of the algorithm and 
the task is to find its asymptotic behavior in the limit of large n ^ 1. The numerator in ® is of 



5 



the order of the largest eigenvalue of dH/dr = Hp — Hb + (1 — '^t)He, which typically scales 
polynomially with n. However, AE'niin can scale down exponentially with n and in such cases the 
required runtime of the quantum adiabatic algorithm to find a solution grows exponentially fast 
with the size of the input. 

One should note that the second term in the r.h.s. of Q is zero at r = and r = 1. Therefore, 
by using different driver Hamiltonians He one can design a family of (possibly random) adiabatic 
evolution paths that start at r = in the same generically chosen initial state and arrive at the 
ground state of Hp at r = 1. In general, different paths will correspond to different minimum 
gaps (yfjnin and one can introduce the distribution of minimum gaps. This distribution can be used 
to compute the fraction of the adiabatic evolution paths / that arrive at the ground state of Hp 
within polynomial time, 

T<cn-", a>0, c = (8) 

For a successfully designed family of paths the fraction / is bounded from below by a polynomial 
in 1/n which leads to the average polynomial complexity of QAA. 



III. BINARY OPTIMIZATION PROBLEM WITH SYMMETRIC COST FUNCTION 

Consider a binary optimization problem defined on a set of n-bit strings z with the cost function 
Ez in the following form: 

n 

Ez = f (Wz) , U!z = ^ Zj. (9) 

This cost is symmetric with respect to the permutation of bits, it depends on a string z only through 
the number of unit bits in the string Wz (the Hamming weight). In this paper we consider the cost 
function Q in the following form which is generalization of the cost introduced in flQ,!^ 

Ez = ^ c{zi^ + Zi^ + Zi^), (10) 

ii<i2<i:i 

C{m) = Podrafi + Pl5m,l + ^2^771,2 + Psf^m.S- 

Here the sum is over all possible 3-bit subsets of the n-bit string z. A subset Zi-^ + Zi^ + Zi^ 
contributes to the total cost a weight factor where A; is a number of units bits in the subset. A set 
of weights [pk] defines an instance of this generalized Hamming Weight Problem (HWP). One 
can formulate a random version of HWP, e.g., by drawing numbers [pk] independently from a 
uniform distribution defined over a certain range. 
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In the limit of large n ^ 1 the cost function dTUt takes the following form: 

3 

E, = PGp{l-^y Gp{q) = Y,(3kq\ (11) 

fc=0 

here I = n/2 and we only keep the terms of the leading order in n. The coefficients Pk in (fTTt are 
linear combinations of 

= I [Pl + (-1)^2] + ^ h + (-1)^3] . (12) 

here = 1 for A; = 0, 1 and = —1 for k = 2,3. 

The function Gp{q) in dTTT) is a third degree polynomial in g, and the form of the function 
depends on the coefficients Pk (pk)- It is easy to show that there is a finite size region in the 
parameter space {pk} where Gp{q) is a non-monotonic function of q that has global and local 
minima on the interval g G (— 1 , 1) . Those minima are separated by a finite barrier with width 5q = 
0{1). The barrier separates strings that have close values of the cost E^^ but are at large Hamming 
distance from each other: they have 0{n) distinct bits. This property can lead to exponentially 
small minimum gaps in QAA due to the onset of low-amplitude quantum tunnelling 141. 



IV. CONSTRUCTION OF THE CONTROL HAMILTONIAN 

A. Representation in terms of operator components of the total spin 

It is natural to consider the control Hamiltonians © for solving the HWPs that are symmet- 
ric with respect to permutation of individual bits Q- In what follows, we use the normalized 
components of the total spin operator S for the system of n individual spins- 1 

1 1 " n 

fij = -Si, Sj = -y a], I = —. (13) 

1=1 

Here Sj are the projections of the total spin operator on the j-th axis (J = x, y, z) and a* are Pauli 
matrices for the i-th spin. For the sake of bookkeeping, in (fT51) and also throughout the paper we 
use "hats" for the spin operators, such as Sj, hj, and some others, in order to distinguish them 
from their corresponding eigenvalues {Sj and rij, respectively, in the above example). 

To obtain the problem Hamiltonian ^ we make use of the obvious connection between the 
values of the Hamming weight function of an n-hii string z and corresponding eigenvalues 
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FIG. 1: Plots of the cost function (fTTl) Gp vs q for different choices of the weights {pk}- Curve 1 corresponds 
to po = 0, pi = 3, p2 = 1; P3 = 1> and the cost function Gp{q) has a global minimum at g = 1, 
corresponding to the string z with the Hamming weight zero, zi = ^2 = • • • = -^n = 0. It also has a local 
minimum at g = — 1 corresponding to the bit string with Hamming weight = Z2 = ■ ■ ■ = Zn = 1. The 
curve 1 yields the particular form of the cost function Gp{q) considered in jj], fl]. Curve 2 corresponds to 
Po = 0.5, pi = 2.5, p2 = —2, p3 = 0.3, it has a global minimum at q = q* inside of the interval (—1, 1). 
This minimum corresponds to approximately bit strings z that all have the same Hamming weight 

Wz = w* = n(l - q*)/2. 

of the spin projection operator 

n^\z) = n^\z), n^ = l-^. (14) 
Then from Eqs. (l2h. (fmi and (fT4ll we obtain 

Hp{n,) = l^Gp{n,). (15) 
We chose the driver Hp in a bit-symmetric form that coincides with © (up to a constant term) 

Hp = P Gp{n,), Gb{x) = 2(1 - x). (16) 
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B. Bit-symmetric drivers He 

It was proposed in [Q] that He can be constructed using some generic ensemble of random 
matrices. The bit-symmetric random drivers for the cost functions of the type (flOb can be con- 
structed as follows 111. One generates an 8 x 8 random Hermitian matrix A with zero diagonal 
elements and non-diagonal elements that are independent random numbers identically distributed 
in a certain interval. Matrix elements of can be enumerated by all possible configurations 

of a 3-bit string {zi.Zj^Zk}- Then He takes the form 

i<j<k zG{0,1}" 

IZijfc) = ® ■ • • ® \Zi)i \Zj)j (g) ■ • ■ (g) \zj,)k \Zn)n- (18) 

Here |z) are computational basis states Eq. ([T]) corresponding to bit-strings z = {zi, . . . , 
and string z-j^ has three of its bits flipped at the positions i,j and k as compared to the string 
z (i.e., Zi = 1 — Zi, etc). Each randomly selected A generates He and therefore a random path 
modification of the QAA. 

From the above discussion, it follows that the matrix of the operator He (flVt is symmetric 
with respect to the bit permutations and therefore it commutes with the operator of a total spin 
S"^ of a system of n spins \. This means that He acts independently in each of the sub-spaces 
corresponding to certain values of the total spin < / < | J22I] . It follows from dTSt and ([T6t that 
the same is true for the total control Hamiltonian © 



Hi 



T 



0, rG(0,l). (19) 



Since in our case the initial state Q is a totally symmetric combination of all states and therefore 
corresponds to the maximal spin / = |, our system always stays in this sub-space during the 
algorithm. Therefore in the analysis of the complexity of QAA one can reduce the 2" x 2" matrix 
of H(t) to the (2n+l) x (2n+l) matrix that only involves the states with different spin projections 
of the maximum total spin / = f ■ Binary strings corresponding to the quantum states from this 
subspace are distinguished from each other by their Hamming weight only. 

In Appendix|Al we show that in the case of real- valued symmetric matrices A and in the large- 
spin limit, the bit-symmetric driver He (fTTt can be presented as a linear combination of 6 operators 
expressed in terms of the large spin operator components n^., acting in the subspace with / = f . 
Using this fact, and also Eqs. dTSt and ([T6b one can write a bit- symmetric control Hamiltonian Q 
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in the following form 
1 



^^H{T) = G{r,hx,h,) = {l-T)GB{nx) + r{l-T)GE{nx,h,) + TGp{n,), (20) 



GEifix, n^) = 7in^ + 72^^ + 73^^ + 74^x^z + l^fixfil + ^ehlh^, (21) 

where {7^} (k = 1, . . . , 6) are independent real coefficients given in Eq. (IA8I) . As we show in 
Appendix |X1 any random realization of the real matrix A can be mapped onto combinations of 
drivers d^TT) by the appropriate choice of the real coefficients 7^. 

We note that He in Eq. ^TH does not have any terms involving Uy operator. The reason for 
that is that we chose matrix elements of A^^zjz^ to be real numbers. Then matrix elements 
of He in z basis are real as well. In this case He can only involve terms with even powers of 
Uy. In (|2ni we have used a conservation of the total spin ( see discussion above) and substituted 
nl = l-nl-nl. 

The form of the total Hamiltonian in (l20b allows us to analyze the minimum gap in QAA with 
different paths © using the WKB analysis of the dynamics of a spin-| in the large spin limit 
(n > 1). 



V. ADIABATIC EVOLUTION OF A LARGE SPIN 

A. WKB approximation for the large spin 

Our analysis in this section is a particular application of the WKB-type approach commonly 
used for the description of quantum spin tunnelling in magnetics 

[Q], S, H, 0. This 

approach is applicable for the large spins (/ >> 1), which is the case of interest for us. 

We choose 2; as a quantization axis and following the standard procedure to obtain the effective 



quasi-classical Hamiltonian in po 



use of the Villain transformation 11231] 



ar coordinates {6, ip} with 9 G [0, vr] and (p G [0, 2n]. We make 



fix = ^/l + e-n^in^ + e) cos (<^) , e = y, (22) 



where azimuthal angle operator (p satisfies the commutation relation 



[(p, hz] = ie. 



(23) 
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In a change of notation we introduce a coordinate q and canonically-conjugate momentum p (cf. 



d 

q = n^, p= -te— = -f, 
aq 

(—1 < g < 1). Expanding (|22t in the large spin limit e ^ 1, we obtain 

(l-q'^Y^^cosp + e- ^^^^ 



2(1 + g) 

Finally, we write the scaled Hamiltonian of the system (l20b in terms of the new variables 

G(n^,n^,r) = n{q,p,T), 

^(?,]5,r) = H(g,p,r) + AH(g,p,r), 



where 



H(g,p,r) = G (^v^l - g2 cosp,q,T 



[q,p] =ie, 



and AH is a small correction 



AH(g,p,r) = €■ 



cosp dG , ^^^2^1 



2(1 + g) (9n, 

(here dG/dn^ has the same arguments as G in (1271 ). 

The stationary Schrodinger equation ^ in the new basis 

7^(g,p, r)^fc(g; r) = Afe(r)^fc(g; r), 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



can be solved in the WKB approximation with the small parameter e ^ 1 playing the role of a 
Plank constant. Then the wave function \l/(g) takes the form 



^(g) = B{q) exp 



iA{q) 



(31) 



where in the leading order in e the action function A{q) satisfies the following Hamilton- Jacobi 
equation 

(32) 



HU,^,H^A. 



dq 

This equation describes a ID auxiliary mechanical system with coordinate g, momentum p, energy 
A, and Hamiltonian function H(g,]9, r). Classical orbits satisfy the Hamiltonian equations 



q{t) = Hp(g(t),p(t),r), p{t) = -H,(g(t),p(t), r) 



(33) 
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where Hg and Hp stand for the partial derivatives of H with respect to q and p, respectively. Sta- 
tionary points of the dynamics (g*, p*) correspond to the elliptic and saddle points of the Hamilto- 
nian function 

Hg(g*,p*,r) = 0, Hp(g*,j9*,r) = 0, (34) 
Elliptic points are minima (or maxima) of H{q, p, r) on the (g, p) plane. They satisfy the condition 

= Hpp{q^,p^,T)Hgg{q^,p^,T) - Hgp{q^,p^,r) > 0. (35) 

where Hpp is understood as a second derivative of H with respect to p, etc. Saddle points corre- 
spond toQl <0 in (l35t . 

In the limit e ^ 1 the adiabatic ground state ^o(g; t) (l30b is localized in the small vicinity of 
the fixed points (g*, p*) corresponding to the global minimum of H(g, p, r) at a given value of r. To 
logarithmic accuracy the WKB- asymptotic (OTT) of the ground-state wave function is determined 
by the mechanical action for the imaginary-time instanton trajectory {q(t),p(t)) emanating from 
the fixed point {q*,p*) 



e 

g(-zoo) = g*, p{-ioo) = p*, g(0) = g. (36) 



^o(g;i") ~ ^o(g;g*,p*,^) oc exp 



dtq{t)p{t) 



too 



Integration in (l36b is along the imaginary axis (—ioo, 0). The instanton trajectory obeys Eq. (l33b 
with the boundary conditions given above and t E {—ioo, 0) corresponding to the line of integra- 
tion in (l36b . The choice of the final instant, t = 0, is arbitrary since the instanton trajectory is 
degenerate with respect to a shift of the time axis. 

We note that the WKB asymptotic ( BHt decays exponentially fast as the coordinate g in (1^ 
moves away from its value at the global minimum g* into the classically inaccessible region. This 
corresponds to the growth of the imaginary part of the action in (l36b . similar to the conventional 
quantum tunnelling in the potential. In the vicinity of (g*,p*) the ground-state wave function 
\l'o(g) takes the form similar to that of harmonic oscillator: 



^o(g) = c X exp 



* f^^ Hgp 2\ rnMlSq"^ 



e 



p* Sq — 6q 



2 H 'p'p j 2 £ 



m* = . , ^, (37) 

\Hpp[q^,p^,T)\ 

here c<j^, > is defined in (B51) . Similarly, the energy spectrum in that region corresponds to the 
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classical elliptic orbits with oscillation frequency 

Afc - H(g=,,p*,r) ~ eil* ^A; + , k = 0,l, (38) 

We note that the frequency depends on r and determines the time-varying instantaneous gap 
between the ground and first exited states, A A = ef2^,(r). 

VI. LOCAL AND GLOBAL BIFURCATIONS DURING THE QAA 

It can be seen from Eq. (l27t that the global minimum of H(g, p, r) will correspond to = ifevr 

(fc = 0, ±1, . . .) provided that the following condition holds for all Ux'- 

_ dG {nx,q^{T),T) 

ujx = r 0, (39) 

where the positive and negative signs of correspond to even and odd values of k, respectively. 
The value of in corresponds to the global minimum of the effective potential [/(g, r) 

U{q) = G{^l-q^, q, r), U{q) - U{q,) > 0. (40) 

Under the above condition the Hamiltonian function of the system near the global minimum (g*, 0) 
exactly corresponds to that of the harmonic oscillator with effective frequency fi* (l35b and mass 
(|37l) 

H(5. P. r) = + '"•"•'-><^-^-'-»' . (41) 

2m* (r) 2 

1 , dG(^l- q^,q^,T 

— = -aA^ ^ ^ 

m* ill = U"{q^). 

In the WKB picture the ground state of the system correspond to the particle performing zero-level 
oscillations near the bottom of the slowly varying potential ?7(g, r). There are two types of the 
bifurcations that can destroy the above adiabatic picture: 



1. Local bifurcation 

Assume that at some instant of time r = tq the effective mass m^ij) goes to infinity. In the 
vicinity of this point the Hamiltonian function (l27t can be approximated as follows: 

H(g, P. r) = ^p' - + ^5q' + dos5q + 0{s'l% 

6q = q- q^To), s = t - Tq, (42) 
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where 



ao 



g*(ro)), &o = 7r-i^\/l-9*(^o), 



, d^U 

^0 = A^' "0 = ^^^^ 
oq^ oqoT 

in the above equations all functions are evaluated at the point {q*{To),p = 0). Equation (l42b 



corresponds to bifurcation point 11241] . It can be seen from (l42b that for r > tq the single global 



minimum of H{q, p, r) splits into the two minima with nonzero momenta 

« ± f !Ml^y'\ ,_,„>o. (44) 



ao 

Due to the symmetry H(g,p, r) = H(g, —p, r) the two global minima with nonzero will stay 
symmetric with respect to the g-axis at later times. 

It follows from (l35b . (14111 ') that the linear oscillation frequency vanishes at the bifurcation point, 
^*(to) = 0, however the energy gap AA(ro) ^ 0. By solving the Schrodinger equation ^ at this 
point in the representation of the momentum p one can find the eigenfunctions ^^(p, tq) and eigen- 
values Afe(ro) corresponding to a ID quantum system moving in a quartic potential (cf. Eq. S42\i ). 
This analysis yields an estimate for the value of the gap, and the characteristic localization range 
5p for k{p, To) 

A A ~ e^/^ 6p ~ e^/l (45) 

The size of the energy barrier in momentum p separating the two global minima in (l42b grows 
with time for r — tq and this leads to a rapid decrease of the energy gap. Sufficiently far from 
the bifurcation point, r — tq ^ e^/^, each of the global minima (g*(r),p^(r)) gives rise to its 
own WKB asymptotic (l36b localized at the minimum. The ground state and the first exited state 
correspond to their symmetric and anti-symmetric combinations, respectively 

^fc(g) = ^ (^(g; q*,pt,r) + (-i)'^(g; g*,p;,r)) , 

A; = 0, 1. (46) 

For T — To ^ e^/'^ the tunnelling splitting of energy levels for the symmetric and antisymmetric 
states determines the value of the gap AA(r) and decreases exponentially fast with r — tq. Away 
from the bifurcation region, r — To = (9(1), the gap scales down exponentially with n (note that 
e = 2/n). 
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As a result of the local bifurcation, the purely adiabatic evolution in QAA collapses. The 
amplitude of staying in the adiabatic ground state for r < tq is nearly equally split between the 
states (l46t with the two lowest eigenvalues. In general, this may reduce the probability of finding 
a system in a ground state at r = 1 by a factor of 2. We note that the control Hamiltonian (l20b 
is at most a cubic polynomial in n^, n^, and therefore the number of local bifurcation events 
during QAA is of the order of one. In the worst case they will cause the reduction of the success 
probability in QAA by a constant factor. 

For a given instance of the cost function (fTTT) defined by the coefficients Pk (or pk) the onset of 
local bifurcations d42b depends on the choice of the driver Hamiltonian He (EB- 

There are a number of ways to select coefficients 7fc's in the driver Hamiltonian (1211) to avoid 
local bifurcations during QAA in a broad range of values of the coefficients pk- For example, to 
completely suppress local bifurcations d42h one can keep in (|2T]) only terms linear in hx and set 



The Hamiltonian function H = H{q,p, r) defines a 3D surface over a 2D plane {q,p) and the 
shape of this surface varies with time r. We consider global bifurcations of this surface where 
the energies of its two minima cross each other at some instant of time r = tq while the distance 
between the minima on the {q,p) plane remains finite at the crossing point. For r > tq the 
minima exchange their roles: global minimum becomes local and vise versa. Before and after the 
intersection in the energy space the two minima are uniquely identified with the ground and first 
exited states of the system's Hamiltonian dTTt . The corresponding wave functions \l/o,i(5') are well 
approximated by their asymptotic expressions(l36b. (l37t . 

The small vicinity of the global bifurcation point can be described within the standard 2-level 
avoiding-crossing picture. There ^o,i are given by symmetric and antisymmetric superpositions 
of the WKB- asymptotic corresponding to intersecting minima. The value of the gap changes with 
time as i/c^(r — tq)^ + AAf^j^ where c is some constant and the minimum gap is determined by 
the overlap of the WKB asymptotic. To logarithmic accuracy it is given by the imaginary part of 
the mechanical action (l36t along the instanton trajectory connecting the two minima 



72 = 73 = 76 = 0. 



(47) 



A. Global bifurcation 




(48) 
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lim g(r) = g^'^, lim p(r) = p^'^ (49) 

T— >±«00 T—^±^00 

Here are coordinates of the two minima; H{ql,pl, tq) = H{q'^,pl, tq), and the instanton 

trajectory obeys the Eqs. (l33l) . The analytical expression for the minimum gap was studied in lO], 
^ for the case He = 0, using a simplified version of the Hamming Weight problem dlOb . Below 
we identify certain geometrical properties of the global bifurcations in the case He = that will 
be used later in the selection of the drivers He for the successful QAA. 



U 



1=. 











q 



FIG. 2: The global bifurcation mechanism: the effective potential profiles U {q, r) vs g for r < tq, t = tq 
and r > To are represented by the curves 1,2, and 3, respectively. 



1. The case He = 

Inthecase7j = (j = 1, ...6), the Hamiltonian has a minimum at p^,(r) = vr/c and the value of 
g*(r) corresponds to the global minimum of the effective potential U (g, r) (BUb . We use Eq. (EUb 
and also the condition U'{q^) = to obtain the following equation for g*(r) 

dq^{T) G'p{q^{T)) 



dr {l-T)U"{qJT).T) 



(50) 



This equation holds until the global bifurcation point at r = tq where g*(T) changes discontinu- 
ously in time (see Fig. 2). At the minimum of the potential U"{q^) > and therefore the direction 
of the motion of g*(r) entirely depends on the direction of the "force", —G'p{q^). At r = the 
potential U{q,0) has a unique minimum at the point q = g*(0) = 0. It is clear that with this 
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initial condition equation (l5Dt can lead to a "wrong" minimum of Gp{q) that lies above the global 
minimum, and such cases will give rise to a global bifurcation. This effect is illustrated in FigH] 
where the two different cost functions correspond to the same direction of motion for g*(r). The 
value of (r) may either smoothly approach the global minimum of the cost (curve 2), or move 
toward a "wrong" local minimum (curve 1), leading to the global bifurcation and exponentially 
small gap in QAA. Adding He to the control Hamiltonian can invert the direction of motion of 
g*(r) toward the global minimum of Gp{q). This can be seen from the fact the Eq.lBOt in presence 
of He possesses the additional term 

(1-r) dGE{^l^^,q.) ^^^^ 
U"{q^,T) % 

(here we drop for sake of brevity the argument r in g*(r)). Clearly, the successful Ge should not 
possess reflection symmetry with respect to n^. Therefore we should only select the terms in (OTl) 
that contain odd powers of n^. Taking into account (l47t we arrive at the following form of the 
driver Hamiltonian 

GE{nx, fiz) = ^iUx riz- (52) 

This driver can remove the potential barrier between the two competing global minima of U (g, r) 
by shifting the original minimum at r = towards the true global minimum of the cost function 
Gp{q) (cf. Fig. [T]). In the classical picture (ETl the driver (15^ corresponds to an external field 
parallel to z-axis which can destroy the tunnelling barrier alongthis direction. The mechanism of 
such tunnelling avoidance is similar to the one considered in yj], where the external field gener- 
ated by the driver (l52b compensates the effective field due to the linear term proportional to the 
coefficient Pi in the problem Hamiltonian (fTTT) . 



B. Bifurcation transition to the tunnelling regime 

In general, one can expect that a complete suppression of the tunnelling barrier at all values 
of T requires a certain magnitude (and sign) of the coefficient 74 depending on the choice of the 
coefficients Pk in the cost function Gp{q). 

The transition to the tunnelling regime can be described as an A3 bifurcation point, illustrated in 
Fig. 2. The effective potential U changes parametrically with r, 74 and {Pk}- Near the bifurcation 
point (tc, 74c, qc) ,the potential has the form U = a 5q^ + h 5'^5q^ + c 5q5T where (5r, 5q, ^7 are 
deviations from the bifurcation point in r, q and 74, respectively. The corresponding conditions for 
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the A3 bifurcation point are: 

dU _ dHl_ _ dHl_ _ ^ ^^^^ 
dq dq^ dq^ 

Taking into account (l40b and (I20b. (|2ni .(l52t. the above equation yields 

= (54, 



74c (1 - Tc 



3 [2-r,(2-/?2)]^ 
(1-r,) (3/3i+/?3) + r,/32/?3 



T,(32 - 2 (1 - 

These equations should be solved for 74c and Tc for the given set of the coefficients The 
bifurcation is avoided when 

I74I > l74c|- (55) 



For example, in the particular case of the HWP (fTOb considered in 111, we have 

A = 1/2, /32 = -3/2, /53 = -7/6, (56) 
r, ^ 0.44, ^ -0.95. 

In this case the example of the driver Hamiltonian He that allows to avoidance of tunnelling in 
QAA was given in 111 where the value of 74 = —8 was used. According to (l56b this value is way 
below the critical value 74c. 



1. Numerical Simulations of the bifurcation boundary 

We performed numerical simulations with the effective potential (l40b checking for the onset 
of tunnelling for all pk E [0; 3], k = 1, ...4. The numeric simulations confirm that the situation 
discussed above is typical for the general HWP, implying that (l52t is the only driver term that can 
be fundamentally responsible for the tunnelling avoidance in a general case, if the coefficient 74 is 
defined appropriately. In particular, one of the two drivers (l52b with 

74 > 7c = 4.9 or 74 < —7c = —4.9 

Pk e [0; 3] A; = 0, 1, 2, 3, (57) 



always suppresses tunnelling in the QAA. 
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FIG. 3: The critical value 7c vs domain size L. 

We also solve the Eqs. (15^ numerically for coefficients taking values on a dense grid of points 
in the cube pk G [0; L] (/c=0,l,2,3)- For each size of the cube L we select the point with largest 
value of 74c denoted below as 7c = 7c (L). The results are presented in Fig. 3. The critical value 
7c is monotonically increasing in L, and the dependence is close to linear for sufficiently large 
L, but it is non-linear in the range < L < 3. It can be inferred from Eq. (l54b that a nonlinear 
dependence of 7c on the scale L is due to the fact that the critical time Tc also depends on L. 

The linear dependence of 7c (L) for large L has a simple intuition. According to (fTTT) and (fT2b . 
the magnitude of Gp{q) is proportional to L. According to Eq. dSOb . the maximal magnitude of 
the coefficient /5i presents a "force" that can possibly move a system into the local minimum at 
small r . From (fT^ . we conclude that |/3i|ujax = '"^^^Pk&[-L,L] Pi = 2L. In the limit of large L, the 
role of the driver Gb in (l54b becomes unimportant. Therefore, the only competing terms are the 
driver Ge and the problem Hamiltonian Gp. The term (BTT) generated by He compensates (BTTl 
the "force" /3i when I74I > l/Silj^ax' ^^'^ therefore in this limit we have 

7c(L)= max 74c^|/3iL,,, = 2L. (58) 

One should note that among the effective potentials generated by choosing different {pk}, there 
are two subsets that can be mapped onto each other by means of the mirror reflection about the 
g-axis, r) — > —U{q,t). We note that the same driver He can not simultaneously suppress 
tunnelling barriers in each of the two mutually symmetric potentials: if the tunnelling barriers are 
not suppressed with 7c, they will be suppressed with —7c, and vice versa. This gives a simple 
intuition for the tunnelling barrier suppression boundary (l55t . 

Finally we conclude, that it is possible to indicate the range of value of I74I such that the driver 
Hamiltonian He = 'J^'hx fiz will play the role of a universal driver that guarantees polynomial 



19 



performance of the QAA for all instances of the generalized Hamming weight problem ( fTOI) provi- 
sory to the mirror-reflection symmetry in the possible choice of the cost functions and the common 
normalization factor L. 



VII. PROBABILITY OF SUCCESS OF THE QAA WITH RANDOM PATHS 

Using the analysis from the previous section one can estimate the probability of success for 
the QAA with random paths proposed in 111 . In that algorithm, the ensemble of random drivers 
He was generated using random 3x3 matrices A^^zjz^ dTTt . It is shown in Appendix lAl that 
for the bit-symmetric optimization problem dTOb the above ensemble is identical to the ensemble 
of independent uniformly-distributed random coefficients 7^ (k=l — 6) that appears in the large- 
spin representation of the driver He (ED- Then for any instance of optimization problem in (flOb 
defined by the set of the coefficients {pk} one should compute the fraction/ of the domain of the 
coefficients {7^} where the following conditions are satisfied: 

(i) . Condition for the nonzero effective mass (l39l) . 

(ii) . The condition (1551 for the complete avoidance of the tunnelling barriers in combination with 

Eq. (l54b for the bifurcation boundary. 

Here we compute the fraction / for the particular instance of the optimization problem (flOb con- 



;ompute 

m 



sidered in |lll|4|,|9|]. In this case Eq. (B^ takes the form 



-2 + T [71 + 272^^ + 87372^ + 74^ 

+75g2 + 276n,g] ^ 0, (59) 



where n,j. = ±a/1 — and provides global minimum of U{q, r) (l40b . The effective mass is 
non-zero if lu^ 7^ 0, and ( 159b yields an estimate on the range of {7^} as 

I72I + iTel < 1 + 1/2 (|7i| + 3 I73I + I74I + I75I) • (60) 

Following 111 we assume that the non-diagonal matrix elements A^^z^z^ distributed in the 
interval [—3, 3]. Making use of dASIl . we obtain 

I72I + bel < 16, |7i| + 3 I73I + I74I + I75I < 50. (61) 
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Therefore, the probability that inequality d^Ut is satisfied is estimated as 1-15^/ (50 x 16) ~ 
0.71875. On the other hand, the values of 74 in dASI I belong to the range, —12 < 74 < 12. Using 
the value of 7c ~ -0.95 given in (l56t we estimate the probability of 74 < —7c to be approximately 
equal to ^ 0.46. Making an approximation that the cases when the effective mass is non-zero are 
statistically independent from the cases when 74 < —7c, we obtain the total probability of success 
as Ptot ~0.46 X 0.71875= 0.334 1/3, which is in qualitative agreement with the numerical 
results of 111. This estimate can be generalized to the case when the matrix elements A^- 
distributed in the interval [— L, L] for sufficiently large L > 3. In this case, the probability that 
(l60b is satisfied remains the same, ~ 0.71875, while the probability that 74 < —7c is estimated as 
(4L — |7c|) /8L. With the assumption of statistical independence, the total probability of success 
is Ptot ~ 0.718 X (4L - 0.95) /8L, and in the limit of large L » 1 we have Ptot ~ 0.359 which 
exceeds slightly the value for L = 3. 



VIII. POLYNOMIAL QAA AND CLASSICAL DYNAMICS OF LARGE SPIN 

In absence of tunnelling, the dynamics of the large spin can be characterized by classical equa- 
tions of motion for the spin projections treated as c-numbers in the form lll8l] 



dS 
~dt 



with 



'uj , S 



(62) 



dH 



dH OH dH 



dS K^^x ^^y 

In coordinate form and in terms of the dimensionless spin projections, this yields 



(63) 



dn^ 



dt 

—— = u^n^ - u^n^, (64) 
dt 

driz 
dt 



x"'yi 



where we took into account that since H does not contain the Sy component, 'uj = ^ = 
{uj^, 0, cj^}. In the case when the "effective magnetic field" 'uj does not explicitly depend on 
time, the system (l62t . (l63l) has two independent integrals of motion 
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S' = Sl + S'y + Sl (65) 

J _ ^ /-J ^\ UJxSx + ^zSz 



^ ^ ^ V^x + ^z 

The first (l65t reflects the conservation of total spin and also holds for an arbitrary time-dependent 
field 'uj = whereas the second integral corresponds to the adiabatic invariant of the system 
(l62b . (l63t . Since in our case Zj = ^ is parametrically time-dependent, the adiabatic invariant 
is conserved approximately for sufficiently slow parametric evolution. Note that the adiabatic 
solutions always play the role of "envelope solutions". This means that on average, the spin 
closely follows the adiabatic solution, but there are fast oscillatory-type motions superimposed 
on the slow adiabatic evolution. Basically, the adiabatic approximation in the classical case is 
applicable when the "slow" motion is much slower than the fast oscillatory motion. This exactly 
corresponds to the adiabatic evolution of the spin system in the quantum case 

Making use of (l65t and taking into account that at the instant r = 0, the total spin was parallel 
to the X-axis, we obtain J/l = = 1, or 

'n = —, (66) 

u 

implying that the total spin is always parallel the effective magnetic field 'to. Therefore, the adia- 
batic evolution of the large spin can be simply described as the situation when the spin follows the 
effective field (on average). 

We note that at this level, there is a direct correspondence between the adiabatic classical so- 
lution and the quasiclassical wave functions of the large spin parallel to "n . From (1^ . it follows 
that this direction can be identified with the effective magnetic field = This justifies the 
"variational" approach introduced in yj, |4|], identifying the variational wave functions with the 
adiabatic ground states along the evolution paths when the total spin is parallel to Vl . Therefore, 
one can observe that in the absence of tunnelling, the general HWP is solved essentially by the 
classical paths of the QAA. 

IX. CONCLUSIONS 

n 

We apply the quantum adiabatic evolution algorithms with different paths [1] to the general- 
ized Hamming Weight Problem that corresponds to the specific case of the random Satisfiability 
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problem defined in dTOt . We show that any random evolution path produced by this algorithm for 
the HWP can be obtained by using 6 specific deterministic basis operators with random weights 
and therefore is parameterized by 6 independent random numbers. Therefore, the approach to 
QAA with different paths can still be reduced to the large spin dynamics for the HWP. We show 
that only one of these "generators" can be a "universal" driver fundamentally responsible for tun- 
nelling suppression for arbitrary HWP and therefore the problem of constructing such a universal 
driver reduces to the definition of its weight 74. Due to the possible reflection symmetry of the cost 
function, any particular case of the general HWP can be solved with one of the two values of the 
weight with I74I > 74c, that is by applying one of the two universal path modifications. We analyze 
the nature of the wave functions along the successful paths and show that it is quasiclassical and 
corresponds to the dynamics of a large classical spin. Therefore, we show that the general HWP is 
solved by completely classical paths of the QAA and present a complete characterization of these 
paths. 

We analyzed in details the types of bifurcations of the effective Hamiltonian function H(q,p) 
that lead to the collapse of the adiabatic evolution. The global bifurcations correspond to the onset 
of tunnelling in QAA and lead to the failure of the algorithm. In contrast, the local bifurcations 
while still corresponding to exponentially small minimum gap only lead to the decrease of the 
probability of success by a factor of 2. Since in a given problem function H{q,p) is a low de- 
gree polynomial in its arguments there are only a few local bifurcations possible. However, the 
phenomenon of local bifurcations may become important for more difficult random optimization 
problems. Assuming the number of such bifurcations M is large the probability of success is 
reduced by a factor of 2^^^. For M that scales up with n that would lead to the failure of the 
algorithm. 
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APPENDIX A: SPIN OPERATOR REPRESENTATION OF MATRIX A 

The random real symmetric 8x8 random matrix A introduced in dTTt . describes the "tran- 
sitions" between each of the 2'^ = 8 states for each clause involving 3 bits lit. This matrix has 
(8 X 8 — 8) /2 = 28 independent matrix elements and can be presented in the form 

A = + ^(2) + ^(2), (Al) 

vv'here A'^^^ A^'^\ A^^^ correspond to the transitions involving one, two and three bits, respectively. 
For each realization, we have 

^^'^ = J E (1 + + s , (A2) 

s,s'=±l 

A^^^ = (a+a+ + a, ) 1 ^ 6, (l + s a^) 

^ s=±l 

+ aa/3 + (T^aJ) - J] 6, (1 + s d^) , 

^ s=±l 

+D {a^a^ + CTj^aa^ag") + E (o-j^ a^a^ + a^a^ a^) , 

where 

Oq,, bss', BajS^, Ca/S-y, D^p-f, Eafi-y (A3) 

are the real coefficients. The indices (a, P, 7) G {1, 2, 3} label the bits, aj, cr^ and cr^ are the 
Pauli sigma- matrices of raising/lowering, s-projection and ^-projection, respectively and s = ±1 
is a spin projection variable. Note that the operator | (1 + s cr^) is a projector onto the spin state s 
for the bit k. Clearly, the number of independent parameters in (|A2ll is 3x4 + 6x2 + 4 = 28, 
where the three terms of the sum correspond to A^^\ A^'^^ and A^'^\ respectively. Note that this 
number of parameters equals the number of independent matrix elements of A estimated above. 
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In the matrix form, the representation (IA2I) yields 



A. 



c 



036++ a2&++ a23&+ 0136+ ai2&+ -B 

a23&+ (22^+- til3^+ C* ai2&- 

CI3&+- ai2&+ -D ai6__|_ ai36_ 

ai2&- ai3&- ai6_- 

036-+ 026-+ 0236- 

023 «2&— 

CLzi>-- 




(A4) 



where the vector ^ of 8 basis states is 



+ + +; + + -; +-+; + + -; - + +; - + -; - - +; 



(A5) 



and the lower left portion of the symmetric matrix Ac is obtained by reflection with respect to the 
diagonal. 

The driver He is obtained by summation over all clauses. In doing this summation, we take 
into account that now the bit indices {i,j,k} E {1,2, ...n} run through all n bits, whereas the 
indices (a, (3, 7) G {1, 2, 3} characterizing the realization of Ac still run through the 3 bits (since 
the same realization of Ac is applied to all triples of bits). The driver He is given by 



HE = H^^ + Hf + Hf, 



(A6) 



where h'^e' 1 correspond to the transitions involving one, two and three bits, analogous 
todATTl. 



H 



(1) 



H 



(2) 



a,i,j,k s,s'=±l 



(A7) 



a,l3,i,j,k 



s=±l 



=±1 



H 



(3) 



+ D [a-a^at + afaja^) + E (^+^"(7+ + a" a/a,") . 

One should note that the second term on the r.h.s. for i/^^ gives a contribution, which is diag- 
onal in representation and therefore leads to the effective "re-definition" of the cost function. 
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Following the logic of [ 1], we disregard such terms. Also, the commutation relations between the 
total spin components give contributions ~ 1 // to the effective potential and can be neglected in 
the large-spin limit. Taking this into account, we obtain from (IA7I) in the large-spin limit 




where Ua = Sa/l is a dimensionless spin projection on a-axis and the coefficients {7^} are given 
by 




(A8) 



In particular, the deterministic driver considered in yj] corresponds to 6_+ = 6+_ = {bg} = 

{aa/s} = B = C = D = E = 0,aa = l and 6++ = -b = -2. It follows from (lA8t that in 

this case, the only non-zero coefficient in (IA8I) is 74 = —8. This corresponds to He = —AnSxS^, 
which is equivalent to He = —2n (S'^S'^ + 6*2 5^.) in the large-spin limit according to the above 
discussion. 



APPENDIX B: BIFURCATION POINT ANALYSIS 



Taking into account only the 74 term in He and expanding up to the 4th order, we obtain the 
conditions U' = U" = U'" = for the ^3 bifurcation point {tc, 7c, x} in the form ll24ll (cf. S53h ) 

3 



Te (2/32 + 6/53X) 

6rc/?3 



-(l-^c) 



2x -X^ + 'JcTc [1 - -X' 



-(1 -rj [2 - 3x2 + 7,r, (-3x)] 
1 - Tc) (6x + 37crc) , 



(Bl) 
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which have to be solved for 7c , Tc and x for the given {[3k}. Solving for x, we obtain condition 
(I5?b in the text. 
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